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Abstract 


We study the chromodynamical gauge symmetry in relation to the internal 
spin structure of the nucleon. We show that 1) even in the helicity eigenstates 
the gauge-dependent spin and orbital angular momentum operators do not 
have gauge-independent matrix element; 2) the evolution equations for the 
gluon spin take very different forms in the Eeynman and axial gauges, but yield 
the same leading behavior in the asymptotic limit; 3) the complete evolution of 
the gauge-dependent orbital angular momenta appears intractable in the light- 
cone gauge. We dehne a new gluon orbital angular momentum distribution 
Lg{x) which is an experimental observable and has a simple scale evolution. 
However, its physical interpretation makes sense only in the light-cone gauge 
just like the gluon helicity distribution Ag{x). 
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The spin structure of the nucleon has been a subject of intense debate for about ten 
years. Much progress has been made in both experimental and theoretical frontiers [|l|. 
However, some of the fundamental theoretical issues remain unsettled, as exemplified by a 
number of recent works in the literature. In particular, in analyzing the spin structure of 
the nucleon, color gauge invariance is still the cause of some confusion. In this paper, we 
intend to explore several important and relevant issues in detail. 

To set the stage, let us briefly recall the forms of angular momentum operator in quantum 
chromodynamics (QCD). In Ref. [Q], a natural gauge-invariant expression is introduced 

<^QCD = Jq + Jg ■, (1) 


where 


Jn= drx 


'll)' —'ll) + 'iIj'x X iDil) 


Jg= d^x X X (E X B) . 


( 2 ) 


The quark contribution Jq contains two terms. The hrst term is obviously the quark spin 
as S = diag(a, a) is the four-dimensional generalization of the familiar Pauli spin matrices. 
The second term is the quark kmetic (or mechanical) orbital angular momentum, in which 
the covariant derivative D = —V -|- igA originates from the quark kinetic momentum [Q . 
We recall that in a gauge theory the kinetic momentum appears to be more physical than 
the dynamical (or canonical) one, the latter corresponding to a partial derivative —zV in 
quantum mechanics [Q]. The gluon contribution to the angular momentum Jg contains the 
well-known Poynting vector, E x B, the momentum density of the radiation held. Some 
recent studies in terms of the above form of angular momentum operators can be found in 
Refs. 

The same QCD angular momentum can be written in an “interaction-independent” form[| 


J = J d^x -'ipjj^'ip + 'ip^xxl—ix/l'ip 

+ E X A + Ei{x X \7)Ai 


(3) 


Because some terms contain explicitly partial derivatives and gauge potentials, the above 
expression is not manifestly gauge invariant. Nonetheless, the physical meaning of each term 
appears to be clear. The hrst term is the quark spin, the second the dynamical (canonical) 
quark orbital angular momentum, the third the gluon “spin”, and the last term the gluon 
“orbital” angular momentum. 


^ Of course, it is not really interaction-free because the color electric field still depends on the 
coupling constant g, 
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According to the above decomposition, one can write down a snm rnle for the nncleon 
spin 0 

i = iAE(Q=) + + Aj(Q^) + l;(q^) . (4) 

Here the matrix elements of the individnal operators are dehned in a nucleon state with 
p^^ = (E,0,0,p) and helicity 1/2, e.g., 

J d^x{E X Ay 

The dependence results from the renormalization of the composite operators. One expects 
that L'g, Ag and L'^ are gauge as well as frame dependent. The purpose of this paper is 
to study how the gauge dependence affects the physical signihcance of the individual terms 
in the above sum rule. In the following we will work in the inhnite momentum frame in 
which the angular momentum operators are dehned from the angular momentum density 
/ In particular, the color electric held is now 

In a recent paper by Chen and Wang p], it was claimed that although the individual 
operators in Eq. (H) are gauge invariant, they have gauge-independent matrix elements in 
the nucleon helicity eigenstates. In other words, Ag, and L'g were said to be gauge 
invariant. The main point in Ref. [§] is that the variations of these operators under a gauge 
transformation may be written as a commutator between the QCD angular momentum 
operator and some other operators, which has vanishing matrix element in a nucleon helicity 
state. If correct, the theorem would have warranted a fresh look at the physical signihcance 
of L'g, Ag, and L'g. 

We hnd that the theorem is contradicted by the following explicit calculation. Consider 
an “on-shell” quark in the state of momentum and helicity 1/2. We calculate the matrix 
element of the gluon spin operator Sg = f d^x(E x A)^ in perturbation theory. Choosing 
the light-cone A~^ = 0, we hnd at one-loop level 




^9 



( 6 ) 


where Cp = {N^ — l)/{2Nc) with W the number of colors, and /i^ are the ultraviolet 
and infrared cutohs, respectively. On the other hand, in the covariant gauge we have 

Aj = . (7) 

A similar discrepancy was found upon calculating the matrix element of the same operator 
in an “on-shell” gluon state. 

The reason for the breakdown of Chen and Wang’s theorem is a subtle one. In a canon¬ 
ically quantized gauge theory, transforming a calculation from one gauge to another is not 
trivial. In particular, the Hilbert space in the covariant gauge contains a nonphysical sector. 
It is not the same as that obtained after quantizing in an axial gauge. For a gauge invariant 
operator O, the gauge invariance of the matrix element means 


(P1|0(1)|P1) = (P2|0(2)|P2), 


( 8 ) 
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where 1 and 2 label the same physical state and operator in different gauges. 0{i) has the 
same functional dependence on the gauge potential although the latter is itself a gauge- 
dependent operator. The transformation from one gauge to another is a “superunitary” 
transformation which takes a state in one Hilbert space to another. The transformation on 
the gauge potential operator is 


A{1Y = A{2)^^ + , (9) 

where is in general a quantum operator, not a c number function. For instance, in going 
from the covariant gauge to the axial gauge, the gauge tranformation is 




( 10 ) 


which is a quantum operator because A~^ is. In Chen and Wang’s proof, it was assumed that 
the gauge transformation parameters are c numbers, which is a strong restriction. [However, 
in certain special circumstances which we will not discuss here, the theorem may be correct 
in leading order perturbation theory.] 

Thus the concept of the gluon spin contribution to the nucleon spin is in general a gauge- 
dependent one. This feature is also reflected in the scale evolution of Ag. In the light-cone 
gauge A~^ = 0, it is well-known that Ag evolves according to the Altarelli-Parisi equation 

i, 


dAg{Q^) 

d\nQ^ 


CXs 




( 11 ) 


where Po = 11 — 2nf/3 with Uf the number of active quark flavors. In the asymptotic limit 
> cxo, the gluon spin grows logarithmically. 


) I axial gauge ^ InQ^, 


( 12 ) 


where the coefficient of proportionality is fixed by nonperturbative physics. 

In the Feynman gauge, the evolution equation becomes much more complicated. In fact, 
the following gauge-variant operators come to mix with the gluon spin 


Oi = - j d^xVA'txAa, (13) 

02 = - J d?x gf^^^A+^^ X . (14) 


[There is no ghost operator here because the ghosts do not carry spin.] Denote the matrix 
elements of the above operators in the nucleon helicity states as oi and 02 . A lengthy 
calculation yields the following evolution equation 


d 

dlnQ2 
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where Ca = Nc. Thus to evolve the gluon spin to a new perturbative scale, one needs not 
only the gluon spin at the starting scale but also the matrix elements of Oi and 02- To hnd 
out the asymptotic behavior as > cx), we diagonalize the upper 3x3 mixing matrix. 
The three eigenvalues are Ai = (11/6)Ca —%/3 = /3o/2, A 2 = (17/24)Ca — ri//3, and hnally 
A 3 = (11/24)C'a — From these, we found out that the leading asymptotic behavior of 
the gluon spin in the Feynman gauge is the same as that in the light-cone gauge, 

)|Feynman gauge ^ ln Q • (16) 


Of course, the coefficients of proportionality in the two gauges are different. 

Given that the gluon spin is a gauge-dependent concept, it is remarkable that its value in 
the light-cone gauge can be extracted from the gluon polarization distribution measurable in 
high-energy scattering. What one extracts in those experiments is of course gauge-invariant 
and is in fact the matrix element of the following gauge-invariant non-local operator [pi]|]. 




r dx— r . 

J —oo 2 J—oo 27r 


(17) 


However, the physical interpretation of this operator is in general not obvious. Interestingly, 
in the light-cone gauge = 0 , the above operator reduces to the gluon spin operator Sg. 
This relationship says nothing about the gauge transformation property of the gluon spin; 
it merely means that the gluon spin in the axial gauge can be obtained from the matrix 
element of a gauge-invariant operator. In other words, the gauge-invariant extension of the 
gluon spin in light-cone gauge can be measured. [This situation is similar in spirit to the 
following example of length in special relativity. The proper length of a pencil is clearly 
frame independent. When we say the length of a house in the frame v = 0.9999c is the 
same as the proper length of the pencil, we are not saying that the length of the house is 
frame-independent. Rather, we are saying that the length of the house in a special frame 
can be known from measuring a frame-independent quantity.] Note that one can easily hnd 
gauge-invariant extensions of the gluon spin in other gauges. But we may not always hnd 
an experimental observable which reduces to the gluon spin in these gauges. As far as the 
nucleon spin structure is concerned, however, the gluon spin in the covariant gauge is as 
interesting as its counterpart in the light-cone gauge. 

Finally, we turn to the orbital angular momentum operators in Eq. (^. The role of 
the orbital angular momentum in parton splitting processes was hrst studied by Ratclihe 
In |]T^, Tang and two of us worked out the leading-logarithmic scale dependence of the 
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orbital angular momenta in the light-cone gauge. 


d (L'g\^ asm ( -ICf f 

dlriQ'^ \L'g J 271 \ 





The hrst term on the right-hand side exhibits the ehects of self-generation of the orbital 
angular momenta. The second represents the generation of orbital angular momenta from 
the quark and gluon spin. The above equation leads to some interesting results about the 
spin structure of the nucleon in the asymptotic limit. As we are going to show below, 
however, the actual operator mixing is more complicated than what is shown in the above 
equations although the result in the aymptotic limit remains intact. 
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We note that in general there is an additional operator which mixes with the quark and 
gluon orbital angular momentum operators, 


AL = J X {—gA)Yip . 


(19) 


Therefore, we proceed to calculate the matrix element of / (Y‘x'ip'^{x x (—^in “on- 
shell” quark-gluon-quark state. At the leading-logarithmic order, it contains the following 
scale-dependent term 

^ \nQ‘^u{x2p) (H- - -ln(a;2/a;i) (xi -F X2 u{xip)el . (20) 

This result indicates that the operator that mixes with L'^ and is in fact more complicated 
than the simple guess AL. The most general form is the following non-local operator 

/ d^xx X f{n ■ d^,n ■ c?^t)t/’^7 4 + h.c. (21) 


where and are derivatives acting on -0 and respectively, and f{x,y) is a function 
X and y and takes different forms at different orders of perturbation theory. Therefore, we 
conclude that to evolve the matrix elements of the gauge-variant orbital angular momentum 
operators is extremely complicated in the light-cone gauge.|^ The same statement applies to 
the orbital angular momentum distributions defined in Refs. |T^, 

The evolution in the Feynman gauge is again different. Here we do not have the problem 
of mixing with infinitely many operators. Apart from the quark and gluon orbital angular 
momentum operators and AL, the ghost held also carries the orbital angular momentum 
L^. Thus, a complete evolution equation will contain at least the mixing of L^, L^ AL 
among themselves and with Ag, AS, oi and a 2 - Because of its limited use, we have not 
calculated the full mixing matrix. However, we did perform a few calculations just to explore 
some of differences. We hnd that the hrst entry in the evolution matrix in Eq. (plSf) changes 
from —{A/YCf in the light-cone gauge to —(7^/3 in the Feynman gauge. The evolution of 
L^ does depend on AL 


151 


dt 


OPg 


(--Cf + -Ca) AL - -CfL' + ... 
V 3 8 J 3 


( 22 ) 


Conversely, the evolution of AL also depends on the other matrix elements 


dAL 

dt 


CXg 


(c. + ic.) 


AL - CrL[ + ... 


(23) 


^Note that the light-cone gauge calculations must be supplemented with some prescriptions for the 
light-cone singularities (additional gauge fixing). In our calculation, we have used a prescription 
such that the regularization is independent of the minus component of the momenta flowing through 
the gluon propagators. In other regularizations, such as the Mandelstam-Leibbrandt prescription, 
the result can be different [l^. Of course, for studying ing truly gauge-invariant quantities, all 
prescriptions are equivalent. 
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These equations would be interesting only if we could find ways to calculate these nonper- 
turbative matrix elements in the Feynman gauge. 

If the evolution of the gauge-dependent orbital angular momentum is complicated, how 
about their experimental measurement? Is it possible, for instance, to have a gauge-invariant 
extension of the quark orbital angular momentum measurable in high-energy scattering 
similar to the gluon spin? A gauge-invariant operator that reduces to the quark orbital 
angular momentum in the light-cone gauge has been discussed recently in Ref. m- We 
note, however, that non-local operators with dependence on spatial coordinates have not 
been seen in factorization of hard forward scattering processes. In particular, inclusive 
deep-inelastic scattering does not depend on these types of operators. 

Given the difficulty of evolving and measuring gauge-dependent orbital angular momenta, 
a question arises naturally as how to incorporate the polarized gluon distribution Ag{x) in 
unravelling the spin structure of the nucleon, particularly since several experiments have 
been proposed to measure Ag{x) in high-energy processes. A satisfactory solution can be 


found by following the approach outlined in Ref. [|^ and taking seriously the suggestion in 

Ref. 0. 

From the off-forward gluon distributions defined from the twist-two gluon operators, one 
can introduce the gluon angular momentum distribution |0 


J,(x) = -x{g(x) + Eg(x)) , 


(24) 


where g{x) is the unpolarized gluon distribution and Eg{x) is the forward limit of an off- 


forward gluon distribution [T^. Jg{x) is gauge invariant, evolves like the twist-two gluon 


distribution, and is accessible experimentally. From this and the gluon helicity distribution 
Ag{x), we can define the gluon orbital angular momentum distribution 


Lg{x) = Jg{x) - Ag{x) . (25) 

Lg{x) is experimentally measurable because Jg{x) and Ag{x) are. The evolution equation 
for Lg{x) is straightforward 

d In ~ + ^)^gn + 

+ + 1 ) - ^Iggi.^)) ^9n + + 1 ) - A'jgg{n)^ , (26) 

where 7 *^ and A'jij are the anomalous dimensions for the spin-independent and spin- 
dependent twist-two operators [Q. However, the catch here is that Lg{x) can be interpreted 
as the gluon orbital angular momentum distribution only in the light-cone gauge. If one 
studies the gluon orbital angular momentum, say in a covariant gauge, Lg{x) would not be 
sufficient. 
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